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We present magnetodielectric measurements in single crystals of the cubic spin- 1/2 compound 
Cu20Se03. A magnetic field-induced electric polarization (P) and a finite magnetocapacitance 
(MC) is observed at the onset of the magnetically ordered state (T c = 59 K). Both P and MC are 
explored in considerable detail as a function of temperature (T), applied field H a , and relative field 
orientations with respect to the crystallographic axes. The magnetodielectric data show a number 
of anomalies which signal magnetic phase transitions, and allow to map out the phase diagram of 
the system in the H a -T plane. Below the 3up-ldown collinear ferrimagnetic phase, we find two 
additional magnetic phases. We demonstrate that these are related to the field-driven evolution of a 
long-period helical phase, which is stabilized by the chiral Dzyalozinskii-Moriya term D M- (V x M) 
that is present in this non-centrosymmetric compound. We also present a phenomenological Landau- 
Ginzburg theory for the MEg effect, which is in excellent agreement with experimental data, and 
shows three novel features: (i) the polarization P has a uniform as well as a long-wavelength spatial 
component that is given by the pitch of the magnetic helices, (ii) the uniform component of P points 
along the vector (H v H z , H z H x , H x H v ), and (iii) its strength is proportional to rj 2 — r/ 2 _/2, where 
7/11 is the longitudinal and r/± is the transverse (and spiraling) component of the magnetic ordering. 
Hence, the field dependence of P provides a clear signature of the evolution of a conical helix under 
a magnetic field. A similar phenomenological theory is discussed for the MC. 

PACS numbers: 76.60. Jx,75.50.Gg,75.25.-j,77.84.Bw 



I. INTRODUCTION 

In recent years there has been an extensive experi- 
mental and theoretical activity in the field of magneto- 
electrics. The magnetoelectric effect (ME) was originally 
predicted by Curie in 1894, 1 and describes the induction 
of electric polarization P by a magnetic field and, vice- 
versa, the induction of magnetization M by an electric 
field E. A phenomenological theory was developed by 
Landau, 2 who considered the invariants in the expansion 
of the free energy up to linear terms in the magnetic 
field. In this approach, symmetry considerations alone 
can fix the form of the coupling between P and M. 3 ' 4 
Based on such symmetry arguments, Dzyaloshinskii pro- 
posed the antiferromagnetic (AFM) compound Cr203 as 
a first candidate for the observation of the ME effect. 5 
Indeed, the electrically induced magnetization (MEs) in 
Cr 2 03 was observed experimentally for the first time by 
Astrov, 6 and soon after, Rado and Folen 7 measured the 
the converse effect, i.e the magnetic-field induced electric 
polarization (ME#). 

Despite the scarcity of compounds with cross-coupled 
magnetic and electric properties, which is related to 
the antagonistic requirements for the existence of both 
orders, 8-10 the research in the field has grown enor- 
mously Since the discovery of ME effects in Cr 2 03, more 
than 100 magnetoelectric compounds have been discov- 
ered or synthesized. 8 The motivation behind this activity 
stems from the novel and fundamental physical phenom- 
ena involved, but also for exciting potential applications, 



including ME sensors and the electric control of mag- 
netic memories without electric currents (and thus Joule 
heating). 8 ' 10 ^ 15 

At the center of interest have been the so-called mul- 
tiferroics, which are materials possessing at least two 
switchable order parameters, such as electric polarization 
P, magnetization M, and strain. 16 ' 17 Multiferroics show 
spontaneous ME effects in addition to those induced by 
external fields. Despite the coexistence of ferroelectricity 
and magnetism, very few multiferroic materials exhibit 
strong coupling between P and M, 10 and in most cases 
the coupling is rather weak. Typical examples are the 
perovskites BiMnC>3 and BiFe03, where the ferroelectric 
transition temperature is much higher than the magnetic 
one. 18 ' 19 

A breakthrough in the field came with the discovery 
of the multiferroics TbMn0 3 2 ° and TbMn 2 5 , 21 where 
the ferroelectricity is directly driven by the spin order. 
The magnetic state in these compounds is either a spi- 
ral (as in TbMn0 3 , 2 ° Ni 3 V 2 8 , 22 and MnW0 4 23 ' 24 ), or 
a collinear configuration (as in Ca3(CoMn) 2 Og 25 ). For 
the understanding of the ME effects in these systems, 
there has been a number of proposed microscopic mech- 
anisms (such as the spin-current model, 3 the exchange- 
striction model, 26 ' 27 and the electric-current cancelation 
model 28 ), in addition to Landau-Ginzburg phenomeno- 
logical theories. 4 ' 16 ' 29-31 

Here we focus on the cubic magnetoelectric insula- 
tor Cu 2 OSeG-3. 32-42 This system crystallizes in the non- 
centrosymmetric space group P2i3, 32 which allows for 
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piezoelectricity and piezomagnetism but not for a sponta- 
neous polarization. The structure is a three-dimensional 
(3D) array of distorted corner-sharing tetrahedra of cop- 
per spin S — 1/2 ions. The unit cell consists of 16 copper 
ions which belong to two crystallographically inequiva- 
lent groups denoted here by Cui and Cu2 (Wyckoff po- 
sitions 4a and 126, respectively). Each tetrahedron com- 
prises three Cu2 and one Cui ion. 

Polycrystalline samples of this compound were stud- 
ied by Bos et al. 33 who reported a transition to a ferri- 
magnetic state at T c = 60 K. Based on neutron powder 
diffraction data, Bos et al. 33 proposed that in this state 
all Cu2 moments prefer to align parallel to each other 
and anti-parallel to the Cui moments. When all Cu spins 
have their full moment, this 3up-ldown state corresponds 
to a 1/2 magnetization plateau, which is realized for H>1 
kOe. 33 The 3up-ldown nature of this plateau was later 
confirmed by high field (14 T) 77 Se Nuclear Magnetic 
Resonance (NMR). 37 

The presence of a finite ME coupling in Cu20Se03 is 
revealed by an anomaly in the dielectric constant, 33,38 at 
the onset of the magnetically ordered state. At the same 
time, high resolution synchrotron x-ray powder diffrac- 
tion data 33 show that the system remains metrically cu- 
bic down to 10 K, with no anomalous change in the lattice 
constant through the magnetic transition. This is also 
supported by infrared 38 and Raman 35 studies in single 
crystals. Furthermore, the NMR study by Belesi et al. 37 
demonstrates that there is not any measurable structural 
distortion even in high magnetic fields. 37 These findings 
strongly suggest that the dielectric constant anomaly is 
not related to any magnetostructural coupling or polar 
structural distortions, 33 but is rather driven by the (pri- 
mary) magnetic order parameter. 

The magnetism of Cu20Se03 below the 1/2 magneti- 
zation plateau is very special. On the theoretical side, the 
reason can be understood from its non-centrosymmetric 
crystal structure that belongs to Laue class T (23). In 
cubic magnets from this class, the magnetic (free) en- 
ergy contains Lifshitz-type invariants, which impair the 
homogeneity of the magnetic ordering and cause a con- 
tinuous, helical twisting of the magnetic order param- 
eter with a fixed sense of rotation, as first predicted 
by Dzyaloshinskii. 43 Microscopically, the origin of these 
couplings relies on the weak (relativistic) Dzyaloshinskii- 
Moriya (DM) antisymmetric exchange that leads to a 
slight twisting of neighboring magnetic moments. 44 ' 45 

The best known cubic chiral helimagnets are the in- 
termetallic compounds and alloys with the "B20" struc- 
ture, which belong also to the Laue class T. 46 In mag- 
netic systems with this crystallographic structure, like 
MnSi, 47 (Fe,Co)Si, 48 , or FeGe, 49 long-period transverse 
flat ferromagnetic helices are observed to form the mag- 
netic ground state. Currently, there is a surge of inter- 
est in these systems because complex topological spin 
textures composed of topologicical solitons called "chi- 
ral Skyrmions" have been predicted to exist in these 
systems under certain additional conditions, 50-52 and 



have been observed in thin films. 53,5 Furthermore, near 
the magnetic ordering transitions, a complex sequence 
of unusual magnetic textures and properties have been 
observed, 55-78 which are usually referred to as the "A- 
phase" . 

Very recently, it has been found that the magneti- 
cally ordered ground-state of Cu20Se03 is in fact a long- 
period helix with a pitch of about 50 nm. The mag- 
netism has also been described as very similar to that 
of the B20 chiral ferromagnets, and the observation of 
field-driven skyrmion textures stabilized in thin film sin- 
gle crystals has been reported by Lorentz-microscopy. 

Theoretically, the chiral long-period modulation is the 
expected behavior of a magnetic system with a sim- 
ple vector order parameter that belongs to one of the 
crystallographic classes that allow for the existence of 
Lifshitz-invariants. The basic spin-structure is dictated 
by the isotropic exchange, and can be well described by 
a collinear ferrimagnctic order with very large exchange 
fields of the order of 50-100 T. 37 ' 79 In fact, detailed micro- 
scopic calculations based on density-functional theory 79 
show that no magnetic frustration should be present 
in Cu2 0Se03, although the magnetic multi-sublattice 
structure appears to allow for geometric frustration. The 
weak relativistic DM interactions are a secondary effect 
that causes a long-period twisting of this primary ferri- 
magnetic order. This is in agreement with the symmetry 
analysis by Gnezdilov et al., 35 who argued that DM inter- 
actions must be important to understand the magnetism 
in Cu2 0Se03. Therefore, the magnetoelectric insulator 
Cu2 0Se03 should be understood as a chiral ferrimag- 
netic helimagnet, and its magnetization process in a field 
is that of the field-driven evolution of a chiral helix. 

Here we report an extensive magnetoelectric study on 
single crystals of CU2OSCO3. Our data demonstrate the 
MEh effect, as well as a finite magnetocapacitance (MC) 
which sets in at the onset of the magnetically ordered 
state (T c = 59 K). We probe these effects in consider- 
able detail as a function of temperature (T) and applied 
field H a . We also study the variation of P and MC with 
respect to the relative orientations of H a with the crystal- 
lographic well as the relative orientation between 
magnetic and ac-electric E fields (for the MC measure- 
ments) . 

Our ME data manifest a number of anomalies which 
signal magnetic phase transitions. From these anomalies 
we map out the magnetic phase diagram of the system in 
the H a -T plane, and demonstrate that there are at least 
two additional magnetic phases below the 1/2-plateau, 
the low-field and the intermediate-field phase. We show 
that in the intermediate phase the field dependence of 
the MEg data is fully consistent with the evolution of 
a chiral conical helix, whose propagation vector is along 
H a . In the low-field phase, the vanishing of the ME re- 
sponse as H — > suggests a multi-domain structure of 
flat transverse helices, where the propagation vectors of 
the different helices are pinned along some preferred axes 
of the system. The transition between the low-field and 
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the intermediate-field phase can then be ascribed to the 
complete alignment of these propagation vectors along 

The above picture of two additional phases at low 
fields seems to be consistent up to 30 K. At higher 
temperatures, further anomalies are observed which sug- 
gest in particular a complex thermal reorientation of the 
anisotropy. We also find a very characteristic double- 
peak anomaly in the MEh effect in a very narrow (2K) 
window close to T c , which might well be related to the 
possible presence of the "A-phase" mentioned above, 
which is typically expected close to T c . 42 

On the theory side, we present a Landau- Ginzburg 
phenomenological theory which captures both qualita- 
tively and quantitative features of the ME# effect, in 
the intermediate and the 1/2-plateau phase. In particu- 
lar, this theory explains the angular dependence of the 
MEh effect, as well as the robust sign change of the po- 
larization P as we go from the intermediate to the 1/2- 
plateau phase. In addition, it provides strong evidence 
that the ME effect in this compound can be attributed 
to an exchange striction mechanism which involves the 
symmetric portion of the anisotropic exchange. Contrary 
to expectations, the influence of the DM coupling in the 
ME effect is heavily diminished by the long-wavelength 
nature of the helimagnetism in CU2OSCO3. 

The theory also predicts that apart from the uniform 
component of the electric polarization (which is what we 
measure in the present experiments) , there is also a spa- 
tially oscillating component (on a mesoscopic scale) that 
is naturally driven by the long-period helimagnetism in 
Cu20Se03. Finally, a similar Landau-Ginzburg theory 
is also given for the MC measurements. 

Our article is organized as follows. In Sec. II we pro- 
vide details on our experiments which were performed at 
EPF-Lausanne. In Sec. Ill we present our magnetization 
measurements. The quasi-static ME and MC measure- 
ments are presented in Sees. V and VI, respectively. The 
theoretical interpretation of the MEh and MC data are 
given in Sees. VIII and IX respectively. Finally, a brief 
summary and discussion of our study is given in Sec. X. 



II. EXPERIMENTAL DETAILS 

High quality single crystals of Cu2 0Se03 were grown 
by the standard chemical vapor phase method. More de- 
tails about the crystal growth can be found in Ref. [37]. 
For our measurements we have used two single crystals of 
Cu20Se03, which we denote in the following by "Crys- 
tal A" and "Crystal B" . Both crystals have a thin rect- 
angular plate shape, with dimensions 1.6x3.6x0.4 mm 3 
(Crystal A) and 1.6x2.3x0.4 mm 3 (Crystal B). The ori- 
entation of the crystal axes with respect to the crystal 
cuts is determined by single-crystal X-ray diffraction. In 
particular, the normal to the widest face of Crystal A is 
[-554], which is about 6° off from the body diagonal [-111] 
axis, while the normal to the widest face of Crystal B is 
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FIG. 1. (Color online) Demonstration of the magnetoelectric 
effect in Cu20Se03 at 15 K. (a) The applied magnetic field 
(right hand axis) is linearly varied with time along the [—554] 
direction, while the electric charge induced on the (—554) 
surface is measured (left hand axis), (b) The induced charge 
from (a) when plotted against the applied field. 



the [100] axis. 

The magnetization measurements were performed in a 
Quantum Design Magnetic Property Measurement Sys- 
tem (MPMS). For the dielectric measurements electrical 
contacts were painted using silver paste. The capaci- 
tance of the crystals was measured using a HP 4284A 
LCR meter. For the temperature and magnetic field de- 
pendence of the capacitance we have employed a MPMS 
and a Janis closed cycle refrigerator, equipped with an 
electromagnet and a magnetic field regulator. The latter 
was also employed for the quasi-static magnetoelectric 
(MEh) measurements. 

For the quasistatic MEh measurements, a slowly vary- 
ing magnetic field is applied to the crystal at a fixed tem- 
perature and the induced electric charge Q is recorded 
with an electrometer. This method is described e.g. by 
Rivera in Refs. [80 and 81]. In a typical experiment, 
we begin by measuring the charge Q at zero field for a 
few seconds. Then, we ramp up the field from to 7.8 
kOe at a constant rate (0.08 T/min to 0.4 T/min) while 
recording the induced charge. The field is then kept at 
7.8 kOe for a few seconds and then the magnetic field is 
decreased to zero, again linearly with time. The same 
procedure can be continued for negative values of the 
field. An example of this field cycling is presented in 
Fig. 1(a), which shows an experiment at 15 K with the 
magnetic field along the [-554] direction and the induced 
charge Q measured on the (-554) surface of Crystal A. 
Figure 1(b) shows the induced charge Q as a function 
of the applied field. A peak, accompanied by hysteresis, 
occurs at 700 Oe. The position of this peak corresponds 
to the field-induced transition observed in the magneti- 
zation process, and will be discussed in detail in Sec. V A 
below. Finally, the quasistatic MEh measurements have 
been performed in the temperature range 10-60 K. 
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FIG. 2. (Color online) (a) The field dependence of the magne- 
tization at 5 K in our single crystals. The external magnetic 
field H a is applied along the [-554] and [100] crystallographic 
directions, which are perpendicular to the widest rectangular 
crystal faces, (b) The first derivatives of the magnetization 
loops. 



III. MAGNETIZATION MEASUREMENTS 

The magnetic properties of Cu20Se03 were first stud- 
ied by Bos et al. 33 in polycrystalline samples. Since 
then a number of magnetization data in single crystals of 
Cr^OSeOs were reported in the literature. 39 ' 40 

Here we would like to map out the H a -T phase diagram 
of Cu 2 OSe03, so that we contrast it to the one obtained 
from dielectric measurements below (Sec. VII). To this 
end, we have taken detailed magnetization measurements 
as a function of T and H a . In addition, we have data for 
two different crystallographic orientations, by applying 
the field H a perpendicular to the widest crystal faces 
of the two samples, namely to the [-554] direction for 
Crystal A and to the [100] direction for Crystal B. 

A representative set of data for these two crystallo- 
graphic orientations is shown in Fig. 2(a), at 5 K. At low 
fields the magnetization varies linearly with H Q , while 
a kink is observed at H c ~ 600 — 800 Oe, which is 
most clearly seen by a peak in the first derivative of 
the magnetization (see Fig. 2(b)). This peak indicates 
a field-induced transition, which was originally reported 
in polycrystalline samples, 33 and was also observed in 
single crystals as well. 39,40 

At higher values of the magnetic field (H a =1.5-2 kOe 
at 5 K), the magnetization reaches a value of 0.53 \ib per 
Cu for both crystallographic directions. This moment is 
consistent with the 3up-ldown ferrimagnetic state that 
was originally proposed by Bos et a/., 33 and recently con- 
firmed by high-field NMR measurements. 37 In agreement 
with previous results, 39 the magnetization reaches this 
1/2-plateau earlier when the field is applied along a cube 
edge, as compared to a space diagonal direction, i.e., 
H [ioo] <h [- 554 1 

plateau plateau' 

The evolution of the magnetization process with tem- 
perature is presented in Figs. 3(a) and (b), for the [-554] 
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FIG. 3. (Color online) M-H curves collected for a number 
of temperatures and with the applied magnetic field H a along 
the [-554] (a) and [100] (b) crystallographic directions. Panels 
(c) and (d) show the derivative dM/dH for three representa- 
tive curves in (a) and (b) respectively. Here we have corrected 
for the demagnetizing field Hd, i.e. H=H a — H^ is the internal 
field. 



and [100] crystallographic orientations, respectively. To 
highlight some of the features we also show in (c) and (d) 
the differential susceptibilities dM/dH. Here, all data are 
plotted against the internal magnetic field, H=H a — H^, 
where the demagnetizing field H^ has been calculated 
following Ref. [82]. 

To map out the H-T phase diagram of Cu20Se03 we 
follow the T-dependence of the two main characteristic 
anomalies in the magnetization process, namely, the po- 
sition H c of the field-induced transition and the onset 
of the 1/2-plateau H p i atcau . Collecting H c and H plat cau 
for both crystallographic directions, we obtain the H-T 
phase diagram shown in Fig. 4. 

We first remark that there are two temperature regimes 
with qualitatively different features, namely above and 
below 30-35 K. In the low-T region we find that 

Hi" 554 ) <Hir, while HgSL < Hg 5 e i, similarly to 
the 5 K data discussed above. These findings suggest 
that the magnetization prefers the body diagonals in the 
low-field phase, but the [100] axes above the file-induced 
transition. 

In the high-T region, we can first observe that 
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FIG. 4. (Color online) The magnetic phase diagram of 
Cu2 0Se03, as follows from the anomalies in the magneti- 
zation measurements shown in Fig. 3. 
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lous T-dependence. The second observation is that, in 
addition to H c and H p i a t C au, there appears a third char- 
acteristic anomaly, denoted by H[3 00 '. This corresponds 
to a second peak in dM/dH which is visible for T>27 
K (Fig. 3(d)). It is noteworthy that in the same T- 
range, the [100] data show a more pronounced hysteresis, 
which is absent from the [-554] data. We also note that 
jjlioo] follows quite closely the T-dependence of He 554 ' 
(see Fig. 4). Altogether, the data clearly reveal different 
anisotropic properties compared to the low-T region, and 
suggest either a coexistence region ("phase separation") 
between H c and H c /, or a separate phase whose nature 
needs to be investigated further. 

The magnetic properties in the low-T region are read- 
ily interpreted by the transformations of a ground-state 
composed of different domains of flat ferrimagnetic he- 
lices. These helices have their propagation directions 
along the easy axes of the cubic system, i.e. the magnetic 
moments are rotating in the plane normal to such an axis. 
Since the single-ion magnetic anisotropy is absent for in 
the present S=l/2 compound, the leading anisotropy in 
this cubic system should be of the exchange type. De- 
pending on the sign of this anisotropic exchange, the easy 
axes are along [111] or [100]. 46 ' 83 The ground-state may 
form a thermodynamical multidomain state, composed 
of these helices, owing to magnetoelastic effects, similar 
to stable domain states in ordinary antiferromagnets. 

First-order magnetization processes transform this 
multidomain state into a single-domain of a conical helix 
that propagates in the direction of the applied field and 
displays a net magnetization. 83 The details of this pro- 
cess depend on the orientation of the field with respect 
to the easy axes and can be hysteretic. As will be shown 
below by considering the possible ME coupling terms in 
the phenomenological theory, for a cubic magnetoeletric 
helimagnet like Cu20Se03, the conical helices display a 
net dielectric polarization, while the effect from the flat 
helices averages out. Therefore, the magnetically driven 



reorientation process may proceed via an intermediate 
domain state that is composed of flat helices without a 
net dielectric polarization and the polarized field-driven 
conical helix. In that case, it is the long-range electro- 
static effect that stabilizes a multidomain state of mag- 
netic helices. 

The magnetization curves suggest that the transfor- 
mation process into the single-conical helix is concluded 
close to the peak of the dM/dH curves, see Fig. 3 (b). For 
higher fields, the single-domain conical helix evolves to- 
ward the collinear ferrimagnetic state, which is achieved 
by a continuous process at i? p i atoau . The T-dependence 
of ifpiateau depicted in Fig. 4, which should trace the 
evolution of the magnetic anisotropy with temperature, 
is clearly anomalous. At low temperatures the effective 
easy axes in the cubic system should be of [100] type 
at high fields. As there exists a net dielectric polariza- 
tion in the conical and collinear phase, this anisotropy 
may be affected by the magnetoelectric couplings and 
may not reflect the purely magnetic anisotropic cou- 
plings which control the behavior of the flat helices at 
low magnetic fields. Our survey of the magnetic phase 
diagram suggests that the magnetic structure is simple 
for fields H > H c as being that of a conical helix that 
is transformed into a collinear state by increasing mag- 
netic field at low temperatures, where also a substantially 
anisotropic magnetic behavior is observed. However, the 
magnetic spin-structure and transformation processes at 
higher temperatures and close to the magnetic ordering 
transition will require further experimental and theoret- 
ical studies. 



IV. DIELECTRIC CONSTANT e r 

The T-dependence of the dielectric constant e r of 
Cu2 0Se03 has been reported by Bos et al. 33 from capac- 
itance measurements on polycrystalline samples, and by 
Miller et al. 38 following an analysis of far- infrared mea- 
surements on single crystals. 

Figure 5 shows our own data from capacitance mea- 
surements in Crystal A with a parallel plate geometry. 
The experiment was performed at zero magnetic field 
with an ac-electric field along the [-554] direction. The 
frequency was 9 kHz and the voltage 1 V. We have re- 
peated the measurement at various frequencies in the 
range 20 Hz-1 MHz but no significant frequency depen- 
dence was observed. The dielectric loss factor (tan J = 
e'r/e'r) was of the order of 5 x 10 -4 in our measurements. 
The T-dependence of e r was also measured for Crystal B 
but the results are identical. 

Our single-crystal data are in good qualitative agree- 
ment with the results reported in Refs. [33] and [38]. 7 
The most important feature is the enhancement of the di- 
electric constant below ~60 K (see Fig. 5), which marks 
the onset of an additional contribution which superim- 
poses the high-T curve. Since this enhancement co- 
incides with the onset of magnetic ordering, it signals 
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FIG. 5. The temperature dependence of the dielectric con- 
stant at zero applied magnetic field. 



the presence of a magnetoelectric coupling mechanism in 
Cu 2 OSe0 3 . 



V. ME H effect 

Here we present direct evidence for the MEg effect in 
Cu 2 OSe03. Namely, that an electric charge can be in- 
duced at the surfaces of the crystals by an applied mag- 
netic field (without an applied electric field). We shall 
also demonstrate that this phenomenon sets in at the 
onset of magnetic ordering. 

In the following, the measured polarization actually 
corresponds to P • e s = Q s /S, where S is the surface on 
which we measure the induced charge Q s , and the unit 
vector e s is vertical to that surface. 



A. ME fl effect: T & H-dependence ("Crystal A") 
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FIG. 6. (Color online) (a) Magnetic field-induced electric po- 
larization in the [-554] direction at various temperatures. The 
magnetic field is also applied along the [-554] orientation. The 
inset shows a magnified view of the induced polarization at 
56 and 57 K. (b) The T-dependence of the induced polar- 
ization at the position of the field-induced transition Pfcf^ 

(open circles), and at the plateau -P^fpiateau ( nue d circles) 
taken from (a). Also included the T-dependence of the sec- 
ond power of the magnetization M 2 at the position of the 
plateau for crystal A (half-filled circles) and crystal B (half- 
filled squares). 



1. H || [-554], e 3 || [-554] 

The MEh effect in CU2OSCO3 was measured in the 
range 15-60 K following the quasi-static method de- 
scribed in Section II. Figure 6(a) shows the magnetic 
field induced electric polarization in the [-554] direction 
for various temperatures. The magnetic field H a was ap- 
plied along the [-554] direction which is perpendicular to 
the widest crystal face (see inset of Fig. 6(b)). 

The data at the lowest temperature (15 K) show a net 
polarization only for H a > 350 Oe (see also Fig. 1). This 
finding is consistent with the low-field picture of multi- 
domain flat helices where the net polarization averages 
out. We also observe that the polarization reaches a peak 
at about 700 Oe, which suggests that the field-induced 
reorientation of the domains is concluded at this point. 
This is also consistent with the fact that the position 
of this peak (point "A" in the following) is very close 



to the field-induced transition in the magnetization mea- 
surements (Fig. 2). 

At higher fields, Pi^f 54 ^ drops quickly with H a and 
saturates to a negative value as soon as the system enters 
the 1/2-plateau (Fig. 2). This change of sign between 
"point A" and the onset of the 1/2-plateau is a very 
robust feature in our measurements, and, as we show 
in Sec. VIII, it is in fact a characteristic fingerprint of 
the evolution of a conical helix under a magnetic field. 

At higher T's, the overall field dependence remains the 
same, but the induced polarization is gradually decreas- 
ing in magnitude. A finite polarization can be measured 
up to 59 K (the anomaly at the field-induced transition 
can be observed up to 57 K) which corresponds to the 
magnetic ordering transition. This immediately tells us 
that the magnetoelectric effect is driven by the (primary) 
magnetic order parameter. 
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FIG. 8. (Color online) (a) The magnetic field-induced polar- 
ization along [-554], measured at 15 K and for various direc- 
tions of the magnetic field with respect to the [-554] axis. At 
6 — the magnetic field is parallel to [-554], while for = 90° 
the magnetic field lies on the (-554) plane, (b) The angu- 
lar dependence of the induced polarization at the plateau, 
at "point A", and at "point B". The solid lines are fits to 
Eq. (14). 



ing order free-energy term responsible for the ME# ef- 
fect scales quadratic with the magnetic order parameters 
(Sec. VIII). 



FIG. 7. (Color online) (a) Magnetic field-induced polariza- 
tion in the [-554] direction at various temperatures. Here 
the magnetic field is perpendicular to [-554]. (b) The T- 
dependence of the polarization data shown in (a) at the field- 
induced transition Pj,^ 5 ^ 4 ' (open squares) , and at the plateau 

Mnd 6 piateau ( nue< i squares) . We also include the corresponding 
data (rescaled by a factor of 1/2) from Fig. 6(b) where H a 
was parallel to [-554] (circles). 



Interestingly, a second anomaly is observed in the 
p[ nd 554 ' ) vs. H curves as we approach T c . This anomaly 
(point A' in the inset of Fig. 6) could be observed only be- 
tween 56-57 K. Despite the fact that the induced charge 
is rather small in this T-range, this feature was repro- 
ducible between different runs. The existence of this sec- 
ond anomaly may signal the onset of a third phase in 
a very narrow T- window close to T c , and may well be 
related to the presence of the so-called "A-phase" . 42 

The T-dependence of the induced polarization at the 

f— 554) 

field-induced transition, PLja , and at the plateau, 
Mnd^pfateau' are P resen ted in Fig. 6(b). Apart from the 
sign difference, we can also observe that Pj^piitcau shows 



2. H a _L [-554] 



-554] 



a much stronger T-dependence than P> 



(-554) 



ind, A 



Further- 
more, we should point out here that the T-dependcncc 
of Pjncfpfatcau sca l es quite well with the second power of 
the magnetization M 2 measured at the position of the 
plateau (see Fig. 6(b)). This happens because the lead- 



We have also studied the T-dependence of P^ 554 ^ 
when H a _L [-554], i.e., on the crystal plate (see inset of 
Fig. 7(a)), and the results are presented in Fig. 7. An im- 
portant difference compared to the data shown in Fig. 6, 
is that here the induced polarization has changed its sign 
and the magnitude is halved. This feature will be ex- 
plained below in Sec. VIII. We should remark here that 
the field-induced transition and the plateau are observed 
at lower magnetic fields in the measurements of Fig. 7(a) 
as compared to the data presented in Fig. 6 (a). This is 
simply due to the lower value of the demagnetizing field 
for the in-plane configuration, as opposed to the out-of- 
plane orientation shown in Fig. 6. 

Beside these differences, both Pfcffi and P^ptLau 
show the same T-dependence with that in Fig. 6(b). This 
is illustrated in Fig. 7(b) where we have included the data 
from Fig. 6(b) (rescaled by a factor of -2). 



B. MEjj effect: angular dependence (T=15 K) 

We have studied the angular dependence of the ME# 
effect for crystals "A" and "B" at 15 K. The ME H effect 
is strongly anisotropic and the direction of the induced 
polarization can be controlled by rotating the applied 
magnetic field. 




H (kOe) * (deg) h (kOe) a (deg) 



FIG. 9. (Color online) (a) The magnetic field-induced polar- 
ization along [100], measured at 15 K and for H a rotating on 
the (100) plane (with H a || [010] at 4> = 0). (b)-(c). The an- 
gular dependence of the induced polarization at the plateau 
and at "point A". The solid lines are fits to Eq. (16). 



FIG. 10. (Color online) (a) The magnetic field-induced po- 
larization along [100], measured at 15 K and for H a rotating 
from [100] (6> = 0) to [03-7] [6 = 90°). (b) The angular de- 
pendence of the induced polarization at the plateau, and (c) 
at "point A" . The solid lines are fits to Eq. 17. 



1. "Crystal A" 

The angular dependence of the ME# effect for the (- 
554) crystal plate is presented in Fig. 8. Here we measure 
the polarization induced along [-554] while we rotate the 
magnetic field from out-of-plane to in-plane, as shown in 
Fig. 8. 

The induced polarization shows an anomaly at the 
field-induced transition ("point A") at low fields. The 
corresponding amplitude is maximum when the magnetic 
field is along [-554] . The angular dependence of the in- 
duced polarization at "point A" and at the 1/2-plateau 
can be both captured by our phenomenological theory of 
Sec. VIII (solid lines in Fig. 8(b)). 

Apart from the anomaly at "point A" , a second 
anomaly (marked as "point B" in Fig. 8) is observed be- 
tween 30° and 60° and between 110° to at least 125°. In 
the region between 60° and 110° this feature cannot be 
observed since the field dependence of the induced polar- 
ization flattens out and merges with the saturated behav- 
ior at the 1/2-plateau. The appearance of two anomalies 
(A and B) must be related to a more complex domain 
reorientation process in this parameter region, which re- 
quires further investigations. 



2. "Crystal B" 

The angular dependence of the ME# effect for "Crystal 
B" was studied in two different configurations. In the first 
configuration (see Fig. 9(a)) we measure the field-induced 
polarization along [100] while we rotate the magnetic field 
in the (100) crystal plane. For <fi=0° the magnetic field 
is along the [010] axis. We see that we can switch the 
direction of the induced electric polarization by rotating 
the magnetic field in the (100) plane. In particular, P^ " 1 
almost vanishes when the magnetic field is parallel to the 



cubic edges [010], [001], and is maximum when it is along 
the cube face diagonals. 

The field dependence of Pjnd " 1 i s similar to the re- 
sults obtained from the (-554) platelet (Figs. 6, 7, and 
8). Some finite polarization is observed below the field- 
induced transition ("point A"), while at higher fields the 

slope of the P vs. H curves changes sign and finally Plnj " 1 
saturates as soon as the system enters the 1/2-plateau. 
The angular dependence of P^d ' & t "point A" and at the 
plateau are displayed in Fig. 9(b) and (c) respectively. 
Both curves follow a sinusoidal dependence of the form 
sin(2</>), which can be captured by our phenomenological 
theory (Sec. VIII). 

In the second configuration, we measure the induced 
polarization along [100] while we rotate the magnetic field 
from out-of-plane to in-plane, as shown in Fig. 10(a). The 

polarization P;^!] '' is almost zero when the field is along 
[100], and rises with the angle 8 between [100] and the 
magnetic field. Figure 10 shows the angular dependence 
of P 1 ^d°' ) at the plateau (b) and at "point A" (c), and 
they both scale as sin 2 6. This dependence will also be 
explained in Sec. VIII. 



VI. MAGNETOCAPACITANCE 

The field dependence of the magnetocapacitance (MC) 
was first studied in polycrystalline samples of Cu2 0Se03 
by Bos et al. 33 Here we present our single crystal data for 
magnetic fields between and 4 T and in the T-range 4.2- 
60 K. As we are going to show, the magnetocapacitance 



MC(H a ) = 



C(H a ) - C(0) 
C(0) 



(1) 



exhibits a rather strong T-dependence and quite 
anisotropic properties. In particular, the magnitude of 
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FIG. 11. (Color online) The magnetocapacitance (MC) of 
Cu20Se03 at 4.2 K for various orientations of the applied 
magnetic field H a in respect to the crystallographic axis. The 
inset shows the MC data for Crystal A plotted together with 
the derivative of the magnetization and the induced polar- 



ization P 



(-554) 
nd 



, all measured with the magnetic field applied 



parallel to [-554]. 



MC depends on the orientation of the magnetic field with 
respect to the crystallographic axes but also on the rel- 
ative orientation between the magnetic H a and the ac- 
electric field E ac . 



A. Representative data with H a || E ac (T=4.2 K) 

Figure 11 shows some representative MC data at 4.2 
K for three different directions of H a with respect to the 
crystallographic axes, and with E ac || H a . Significant 
MC is observed which is dominated by a peak occurring 
slightly before the system enters the 1/2-plateau (see in- 
set of Fig. 11). The amplitude of the MC at the position 
of the peak is maximum when H a is along [-554], and 
is practically zero when it is directed along [100]. At 
higher magnetic fields where the system is already in the 
1/2-plateau, the MC saturates to small negative values in 
all directions. In the inset of Fig. 11, apart from the MC 
data for Crystal A we have also included the polarization 
data P;^ 55 ^, as well as dM/dH (from Fig. 2(b)) for the 
same crystal, to demonstrate that the anomalies found 
in the magnetization measurements track the ones ob- 
served in the dielectric measurements, which highlights 
the strong ME coupling in Cu20Se03. 



B. T & H-dependence of MC 

1. H a || [-554], E ac || [-554] 

The MC was first studied for the direction exhibiting 
the maximum effect, namely the [-554] crystal axis. Fig- 
ure 12(a) shows the MC data collected with both H a and 
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FIG. 12. (Color online) (a) The magnetic field dependence of 
the MC at various T's. Both magnetic and ac-electric fields 
were applied along [-554]. The inset shows a magnified view 
of the MC at 4.2 K at low fields, where a small shoulder 
("point A") accompanied with some hysteresis is observed, 
(b) The magnitude of the MC at three characteristic values 
of the magnetic field, i.e., at points "A" and "B", and at the 
onset field of the 1/2-plateau (rescaled by a factor of -6.6). 



E ac along [-554]. By increasing T, the MC at the main 
peak ("point B") is smoothly decreasing in magnitude 
and the position of the peak shifts to lower values of H a . 
Apart from the main peak, a small shoulder ( "point A" ) 
with some hysteresis is observed at 700 Oe, which can 
also be seen in the inset of Fig. 12(a). The position of 
this shoulder is the same with the field-induced anoma- 
lies observed in the ME# and the magnetization mea- 
surements (see inset of Fig. 11). The anomaly at "point 
A" is smeared out (and thus it is not discernible) between 
20-30 K, but it can be clearly seen again above 30 K as a 
small dip and with opposite (negative) MC sign. We re- 
call here that such a different qualitative behavior above 
and below the region 20-30 K was also observed in the 
[100] magnetization measurements (Fig. 4), and thus the 
two effects must be correlated. Finally, the MC at the 
plateau is generally quite small and negative, but picks 
up a small positive value above 55 K, which however does 
not saturate up to 4 T. 

Some of the above features can also be seen in 
Fig. 12(b), which shows the MC at "point A", at "point 
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FIG. 13. (Color online) a) The magnetic field dependence of 
MC at various temperatures. The ac-electric field was ap- 
plied along the [-554] crystallographic orientation, while the 
magnetic field is perpendicular to the electric field, (b) The 
T-dependence of the magnitude of the MC at saturation and 
at the point "B", taken from (a). 

B", and at the onset of the 1/2-plateau (rescaled by a 
factor of -6.6). It is noteworthy that the T-dependcncc 
of the MC at "point B" is very similar to the one at the 
position of the plateau. 



2. H a _L [-554], E ac || [-554] 

The MC was also measured with H a in the (-554) 
plane. Here the ac-electric field E ac is kept in the [- 
554] direction, as shown in Fig. 13. As in Fig. 12, the 
MC is again dominated by a peak ("point B" in Fig. 13) 
slightly before the onset of the 1/2-plateau, where a neg- 
ative MC is observed. At low fields (250 Oe) and for 
T < 23 K, a step-like behavior is observed in the MC 
data ("point A"). The field position of this step coincides 
with the field-induced transition observed in our magne- 
tization measurements (not shown here). This low-field 
step-like shape of the MC curves is smeared out for T> 
23 K. 

The T-dependence of the MC at "point B" and at the 
position of the plateau is presented in Fig. 13(b). The 
overall dependence is very similar to the data shown in 




H a (kOe) * (cleg) 

FIG. 14. (Color online) The angular dependence of the MC, 
measured at 15 K with the ac-electric field along [100] and the 
magnetic field rotating on the (100) plane (at <j) — 0, H a || 
[010]. (b) The angular dependence of the magnetocapacitance 
at "Point B" is given (b). 



Fig. 12. However, the MC at "point B" is lower for the 
in-plane configuration, while larger (negative) MC values 
are now observed at the plateau. 



C. Angular dependence of MC (T=15 K) 

We have studied the angular dependence of the MC at 
15 K for "Crystal B" in two different configurations. In 
the first configuration shown in Fig. 14, E ac is along [100] 
while H a rotates in the (100) plane. Both the amplitude 
and the shape of the MC curves change drastically upon 
changing the orientation of H a . As in the ME# experi- 
ments (Fig. 9), the maximum effect is observed when H ( , 
is parallel to the cube face diagonals. On the other hand, 
the minimum MC value is observed when H a is along the 
cube edges. The MC at the position of the maximum 
("point B") is shown in Fig. 14(b) as a function of the 
angle <j> between H a and the [010] direction. A sinusoidal 
angular dependence of the form ~ sin 2 (20) is observed 
which suggests a free-energy term of the form ~ P 2 M A 
which is quartic in the magnetic order parameter. As we 
explain in Sec. IX, such a term arises naturally in the an- 
alytical expansion of the free-energy in powers of the ME 
coupling mechanism that is related to the ME# effect. 

In this second configuration, E ac is fixed along [100] 
while H a makes an out-of-plane rotation as shown in 
Fig. 15. The MC follows a sinusoidal angular depen- 
dence of the form ~ sin 2 6, both at "point B" and at the 
1/2-plateau. As we discuss in Sec. IX, this dependence 
suggests the presence of the more conventional ~ P 2 M 2 
free-energy term in the free energy. Finally, it is quite 
interesting to contrast the large (and positive) MC mea- 
sured e.g. at "point B" when H a _L E ac (first config- 
uration) to the vanishing MC when H a || E ac (second 
configuration) . 
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FIG. 15. (Color online) (a) The angular dependence of the 
MC measured at 15 K, with the ac-electric field along [100], 
and H a rotating from [100] (6 = 0) to [03-7] (0 = 90°). (b) 
The angular dependence of the MC at the plateau and at 
point "A", as taken from (a). 



VII. PHASE DIAGRAM FROM 
MAGNETIZATION & DIELECTRIC 
MEASUREMENTS 

From the position of the characteristic anomalies ob- 
served in the ME H (Figs. 6 and 7) and MC (Figs. 12 and 
13) measurements we have deduced the T-H a phase dia- 
gram (Fig. 16) of Cu2 0Se03 for both in and out-of-the 
(-554) plane orientations of the magnetic field. For the 
out-of-plane configuration we have also included the data 
points from our magnetization measurements for compar- 
ison. The field position at which the low field-induced 
transition and the plateau are probed from the dielectric 
measurements fits quite well with the corresponding po- 
sitions from the magnetization measurements. The main 
differences between Fig. 16(b) and Fig. 16(a) are due to 
the different demagnetizing field for the in-plane and the 
out-of-plane configuration. 

Altogether, the dielectric measurements are consistent 
with the physical picture we described in Sec. Ill above, 
namely that we have at least two additional phases be- 
low the 1/2-plateau, the low-field (H < H c ) and the 
intermediate-field (H c < H < Hpi ate au) phase. The low- 
field phase is a multi-domain chiral phase whose ME re- 
sponse defines two sub- regions (not specified in Fig. 16): 
The first shows a vanishing ME# effect which is con- 
sistent with the picture of multi-domain averaging of 
flat helices, while the second region marks a complex 
field-induced reorientation process which concludes at 
the "point A" peak anomalies. 

The intermediate-field phase is related to the evolution 
of a single-domain conical helix under a magnetic field, 
with the propagation vector along the field. In fact, the 
phenomenological theory presented in Sees. VIII and IX 
below, shall provide very strong arguments in favor of 
this picture. 

In addition to the above phases, there are indications 
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FIG. 16. (Color online) (a) T-dependence of the field position 
of the characteristic points "A" , "A " , "B" , and the onset of 
the 1/2-plateau, as observed in the field-induced polarization 
(Fig. 6(a)) and MC measurements (Fig. 12 (a)), (b) Same as 
in (a), but the results are now from the polarization data of 
Fig. 7(a) and the MC data of Fig. 13(a). 



in our data for at least one more phase. For example, 
the second anomaly (point A' in Fig. 16(a)) observed in 
the small T-window (2K) close to T c might be related to 
the presence of the so-called "A-phase" that is found in 
similar compounds. 55-78 



VIII. LANDAU-GINZBURG THEORY FOR THE 
ME H EFFECT 

A. Conical helix with k || e^. 

Here we develop a phenomenological Landau- Ginzburg 
approach to account for the ME# effect. We shall focus 
on the intermediate- field and the 1/2-plateau phases. In 
this region, the magnetic state can be described by a 
single-domain conical helix with propagation vector k 
along the field direction . We shall keep track of both 
Cul and Cu2 magnetic order parameters by writing 

m (r t )=rie h +^{^e^+h.c.) (2) 

with i — 1,2. Here ip = (ei — ie 2 )/2, and {e 1 ,e 2 ,e; l } 
form a triad with ei x e 2 = e/j. Since the Cul and Cu2 
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moments are antiparallel, T[{ = —f]\ = m and similarly 
-rji = r)±. The long- wavelength transverse twist- 



ing between two nearby moments i and j is captured by 
the relative phase factor e l]si ' Srij (where <5r.y = — rj), 
with k-^r^j -c 1. 

We continue by recognizing that the magnetic and elec- 
tric polarizations of the system form three-dimensional 
(3D) representations of the present cubic group P2i3, 
but the two quantities transform differently under time- 
reversal. Table I summarizes the transformation proper- 
ties under the three point-group generators of P2i3 plus 
the time-reversal invariance. It is also useful for what fol- 
lows to define the following composite 3D representations 
ofP2i3 



{A,B}+ = (A y B z , A z B x , A x B y ) 

{A,B}_ = {B, A}+ = (A z B y , A x B z , A y B :c ) 



(3) 
(4) 



given two vectors A and B. 

Now, to leading order we may look at the invariant lo- 
cal terms that are linear in P and quadratic in the mag- 
netic order parameters. In particular, the theory may 
contain any pair among rj 1 and rj 2 . For any such pair 
(i and j, in the following), there are only two invariants 
that do not involve spatial derivatives (the latter are de- 
scribed in App. XI A): 

rih z P x + V z r, x P y + rjZrjtP* = R, Vj}+ • P (5) 
v!v V jP x + V*V Z jP y + V!V*P Z = {ri 3 , »?<}+ • P (6) 

where {'H i ,r]j}± are two composite 3D polar representa- 
tions of P2!3, as defined above. So to lowest order, the 
free energy reads 



$({ J? },P) = <I>o({r7}) + ^-P-E 



€i{Vi,Vi}++&{Vi,rii}-) ■ (7) 



Here $o stands for the magnetic portion of the energy 
in the absence of the ME coupling (i.e. it contains the 
strong exchange energy plus the chiral DM term and 
other types of anisotropies), \e is the electric suscep- 
tibility (scalar for cubic groups) in the absence of ME 
coupling, and £i : 2 are the ME coupling constants. The 
physical origin of the last term of Eq. (7) becomes more 
transparent by rewriting it as 



1, 



£+ ({Vi,Vj}++{Vi,Vj}-) +t-Vi x Vj 



(8) 



It is now evident that the ME# effect will be dominated 
by the symmetric exchange component, since nearby 
spins are almost collinear in the long-wavelength coni- 
cal state. Thus, despite the fact that the DM coupling is 
the primary source of hclimagnetism, its influence on the 
MEh effect is heavily diminished by its long-wavelength 
nature. 

To make these points more explicit and obtain the final 
predictions for the angular and field dependence of the 
MEff effect, we proceed by replacing Eq. (2) in Eq. (9) 
One immediate prediction is that the helimagnetism in- 
duces both a uniform as well as a spatially oscillating 
polarization component which is also of mesoscopic na- 
ture. In the following we shall restrict ourselves to the 
uniform component, P unl , which is the one we measure 
in our experiments. Using the relation A{xp,xp*}± — 
— {eh,e/j} + ± ie/j, and expanding e lU ' Srij ~ 1 + k- Srij, 
we find: 



-r)j_ (k • Srij) e h 



(10) 



where cr^ is the sign of r/f rjj . The second term is the con- 
tribution from the DM coupling and can be disregarded 
since k • Sr^ < 1, as discussed above. Similar contri- 
butions that scale linearly with k arise also from the in- 
homogeneous Lifshitz-type invariants in the free energy 
(App. XI A). Altogether we may write 



puni 

ijXe 



1 



{e h ,e h }^ 



(11) 



which is our central expression. As we are going to show 
below, the angular dependence found in the experiments 
is fully captured by the projection of the vector {e^-, &h}+ 
along e s . Similarly, the observed field dependence can 
be fully captured (in the field range where the magnetic 
state is described by Eq. (2)) by the quantity m 



1-T7 2 



This quantity goes from —1/2 for a completely flat helix 
(7711 = 0) to +1 at the 1/2-plateau (where t]± — 0), and 
thus provides a physically transparent explanation for 
the sign difference between the measured polarization at 
"point A" and that at the 1/2-plateau. Of course, the 
"point A" does not correspond to a completely flat helix 
(this would be the case at zero-field) and this is why 
the relative magnitude | -Pa |/|Ppiateau| is actually smaller 
than 1/2. 



where £± = £1 ± £ 2 ■ The first term inside the brackets in- 
volves the symmetric exchange anisotropy, while the sec- 
ond involves the antisymmetric DM interaction. These 
two anisotropy types show drastically different contribu- 
tions to the MEh effect. Indeed, minimizing $ with re- 
spect to P gives (for E = 0): 

P/Xe = % ({r h ,V j }++{Vi,ri j }-) + ^rVi x Vj -(9) 



B. Flat helix with k S/f e h 

Here we comment what happens for a flat helix whose 
propagation vector is not along e^, but along some gen- 
eral direction e&. This analysis explains why the total 
polarization averages out in the low-field region of the 
multi-domain helical phase. 
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TABLE I. Symmetry transformations under the three point group generators of P2i3 and the time reversal T. Here P is the 
electric polarization, while rj stands for any of the magnetic order parameters of the system. 





{C 2 *|(iO,i)} 


{C2 V m |,|)} 


{C 3 \0} 


T 




(-rix,-Vv,+Vz) 










(-P x ,-P y ,+P z ) 


(-P X ,+Py,~P Z ) 


{Pz ; Pc , Py) 


(Pc, Py, Pz) 


(d x , d y , d z ) 


(-d x ,-d y ,+d z ) 


(-d x ,+d v ,-dg) 


(9 2 , d x ,d y ) 


{d x ,dy,d z ) 



Given that ei x &i is now equal to e^, one may easily 
show that the above formulas remain valid if we replace 
by efc everywhere. In particular, the net polarization 
is given by 



oc 



{e fc ,e fe }. 



(12) 



k-d> 



omains 



Now, the magnetization process close to the 1/2-plateau 
(Fig. 2(a)) suggests the (100) as the easy axes, but the 
low-field data cannot confirm this (in fact they seem 
to suggest the (111) as the easy axes but this is most 
likely related to the complex low-field domain reorienta- 
tion process, as discussed above). In any case, one can 
show that the net polarization vanishes for both (100) 
(111) cases. In the first case, e& || (100), all three do- 
mains give a vanishing contribution, while in the second, 
efc || (111), each of the four domains gives a finite contri- 
bution (for a flat helix r)?< — rjj_/2 = —1/2) but the net 
polarization from all domains averages out. 



C. MEh effect: Comparison with exp. data 

1. Angular dependence 

We now proceed to a direct comparison to the exper- 
imental MEjj data, beggining with the angular depen- 
dence. In the following we shall use the scaled quantity 



P 



P 



We begin with crystal A for which e s = 
For the data shown in Fig. 8(b), = cos( 



where e' A 
we find 

p(-554) = 



(13) 



e s +sin# e' A , 



5 
451 



= (-41,-25,-20). Applying Eq. (10) 



55 



The solid line shown in Fig 
sion. Parenthetically, we 
proportional to sin(2#) in 
deviation (6°) of e s from [- 
tionality constant in front 
than the first two terms of 
along [-111] we would get 



150cos(26») - 3V41sin(26») 
(14) 

8(b) is a fit using this expres- 
should remark that the term 
Eq. (14) arises from the small 
-111]. This is why the propor- 
of this term is much smaller 
Eq. (14). Would e s be exactly 



p(-m) = l + 3cos(2fl) 
4^3 



(15) 



These expressions explain in particular why the polariza- 
tion data with H a || e s (9 = 0) are larger by a factor of 
-2 from the data corresponding to H a _L e s (6 — 90°), 
see Fig. 7. 

Let us now turn to Crystal B for which e s = [100]. For 
the data shown in Fig. 9, the magnetic field rotates in 
the plane defined by = cos0 [010] + sin^ [001], and 
Eq. (10) gives 



P( 10 °) = X - sin(2^) 



(16) 



which is also in perfect agreement with the data. 

Similarly, for the data shown in Fig. 10, the magnetic 
field rotates in the plane defined by = cos0 [100] + 
sin# e^, where e' B — -^=(0,3,-7), and Eq. (10) gives 
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p(ioo) = sm - 

58 



(17) 



which is again in very good agreement with the data, 
apart from a very small shift which may originate either 
from a tiny misalignment of e s away from [100], or from 
additional minor contributions to the ME free energy. 



2. Field dependence 

Let us now turn to the field-dependence of the mea- 
sured polarization and compare it with Eq. (10). To this 
end we shall make use of some well known and rather 
simple theoretical expressions for the evolution of ry" and 
rj 1 - as a function of applied field, in the conical phase. 
Following e.g. Ref. 84, we write 



V\\ = .9 



H 



plateau 



v± = a 



l 



H 2 /H? 



plateau 



(18) 



where g is a T-dependent prefactor which vanishes at T c , 
and -f/piateau is the onset field of the 1/2-plateau. 

For concreteness we consider the data shown in Fig. 9. 
According to Eq. (16), 



P( 10 °) _ g 2 / H 2 
aijXe£,+ 4 V H 2 



plateau 



l) sin(20) , (19) 



where we note that the polarization scales quadratically 
with H. From our data we have access to both param- 
eters g and -ffpiatcau (the latter is influenced by the de- 
magnetizing field), and so we can make a direct compar- 
ison between theory and experiment. This comparison 
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H (kOe) 

FIG. 17. (Color online) Comparison between theory (solid 
lines, see text) and experimental data (symbols) for the field 
dependence of the electric polarization along [100] for some 
representative angles 8 (data reproduced from Fig. 9(a)). The 
shaded area shows the low-field phase (which further consists 
of two sub-regions). The first arrow designates the onset of 
the single-domain conical phase (peak at "point A") for 9 = 
45°. This phase evolves continuously toward the 3up-ldown 
ferrimagnetic phase (1/2-plateau). 



is shown in Fig. 17 for a few representative angles 9, 
and is remarkably successful. The deviation around the 
1/2-plateau can be ascribed to thermal fluctuations (the 
experimental data are taken at 15 K while the above ex- 
pression is valid only at T=0). 

Altogether, this gives further confidence that the ME 
coupling mechanism of Eq. (7) above captures the exper- 
imental findings both qualitatively and quantitatively. 



IX. LANDAU-GINZBURG THEORY FOR THE 
MAGNETOCAPACITANCE 

A. Theory 

To capture the physics of the MC, we must include 
higher order terms in the ME portion of the free en- 
ergy, and in particular we must look at quadratic terms 
in the polarization. To satisfy time-reversal invariancc 
this term must scale as (j]) p PP, were p is an even in- 
teger. Now the strong ~ sin 2 (20) angular dependence 
of Fig. 14(b) suggests that we should take p = 4, while 
p = 2 may also be present as suggested by the angular 
dependence of Fig. 15(b). 

The simplest possible term with p = 4, that is also 
physically motivated, is the square of the symmetric com- 
bination of the MEh invariant of Eq. (7), namely 

U = ^[({r, i ,r, j } + + {r li ,T ]j }-)-P] 2 = ({r, 1 , Vl } + -P)\ 

(20) 

where in the second equality we disregarded terms re- 



lated to the small phase difference k • 5ry, as discussed 
above. Physically, Uq stands for the second-order term 
in the analytical expansion of the free energy in the ME 
coupling mechanism discussed in Sec. VIII. Including this 
term gives 



S$ = $' - $ = -( U , 



(21) 



where Co is the coupling constant, and $ is given in 
Eq. (7) above. 

The polarization is again obtained by minimizing with 
respect to P. However, in the experiment one uses an 
alternating ac-electric field and measures the alternating 
portion of the charge on the capacitance plates. So from 
the derivative of $' with respect to P we can disregard 
the terms that are related to the MEjj effect. This gives 



T7 n 



pa 

ac 



Co 



dU 



With 



Xe dP a 
E/3 X~p p L we get 



a/3 



d 2 U 

QpaQpP 



(22) 



(23) 



Now, given that the magnetocapacitance is a very small 
quantity (MC <C 1%), we may write XeX 1 — l~$x/Xe, 
and thus 5x/Xe — — Xe^X -1 - If we are measuring the 
charge Q s on a surface S (which is perpendicular to the 
applied ac-electric field), then 



MC = 



5Xs 

Xe 



a/9 



P. 



d 2 Up 
dP a dPP ' 



Taking derivatives we find 
d 2 U 



dP a dPP 



(24) 



(25) 



For example, for a = /3 = x (which is directly relevant to 
most of the MC data we show here), 



d 2 u a 
<)./->- 



2(%V) 



Taking the uniform portion of (rjiTjf) for the conical 
helix gives 



where 



MC 

/i(H„) 
/a(H„) 
/s(H„) 



h KelY + h (e x h Y+h 



_L\4 



_L\4 



±\2 



2 Wii ) 



(26) 

(27) 
(28) 
(29) 



B. Angular dependence of MC: Comparison with 
exp. data 

We now consider the data shown in Fig. 14(b) for the 
angular dependence of the MC at "point B", when the 
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field rotates in the (100) crystal plane. Replacing = 
cos0[O1O] +sin0[OOl] in Eq. (26) yields 

^^Asin 2 (20) + / 3 , (30) 

which is exactly the angular dependence observed in 
Fig. 14(b). This provides confidence that the invariant 
IAq must be present in <5$. 

Next, we consider the angular dependence shown in 
Fig. 15(b). Replacing e h = cos 0[1OO] + sin</>[03 - 7}/V58 
in Eq. (26) now gives 

MP 01 2 

^ = ^/i S m 4 + / 2 cos 2 + / 3 . (31) 

At the 1/2-plateau, fa = 1, fa = fa = (since r]^ = 0), 
and so this expression predicts that MC oc sin 4 0, which 
is clearly in disagreement with the ~ sin 2 behavior of 
the data. So we may conclude that there is at least one 
more invariant giving a ~ sin 2 contribution. In particu- 
lar, this term must scale quadratically with the magnetic 
order parameter (p — 2) to account for the angular de- 
pendence. A list of such invariant terms is provided and 
discussed in App. XI B. 

X. SUMMARY AND DISCUSSION 



effect which might be due to the possible presence of the 
"A-phase" . 

We have also developed a Landau-Ginzburg phe- 
nomenological theory which accounts for the ME# effect 
in the intermediate-field and the 1/2-plateau phase, and 
explains the angular, temperature and field dependence 
of this effect in a simple and physically transparent way. 
Contrary to expectations, the influence of the DM cou- 
pling in the ME effect is heavily diminished due to the 
long-wavelength nature of the twist. Instead, the ME ef- 
fect is driven by an exchange striction mechanism involv- 
ing the symmetric portion of the exchange anisotropy, 
with three central features: (i) the uniform component 
of P points along the vector (H v H z , H z H x , H x H v ); (ii) 
its strength is proportional to i]^ — ?7j_/2, where jyy is 
the longitudinal and r)± is the transverse (and spiraling) 
component of the magnetic ordering. Hence, the field 
dependence of P provides a clear signature of the evolu- 
tion of a conical helix under a magnetic field; and (iii) 
apart from its uniform component, the polarization P has 
also a long-wavelength spatial component that is given 
by the pitch of the magnetic helices. This effectively 
mesoscale antiferroelectric structure could in principle 
be probed by experiments with high spatial resolution, 
which may further clarify the ME coupling mechanism 
in this Cu20Se03 at a more microscopic level. 



We have observed the ME effect in single crystals of 
the cubic compound Cu 2 OSe03, and we have studied its 
angular and temperature dependence. We have demon- 
strated that the electric polarization can be controlled by 
an external magnetic field. We have also investigated the 
temperature, field and angular dependence of the mag- 
netocapacitance. Both the magnetic field-induced polar- 
ization P and the magneto capacitance MC set in at the 
magnetic ordering temperature and show a number of 
anomalous features which we ascribe to magnetic phase 
transitions. From the positions of these anomalies we 
have mapped out the phase diagram of the system in the 
H a -T plane. In the field region below the 3up-ldown 
collinear ferrimagnetic phase and for T <30 K we find 
two additional magnetic phases, whose origin is related 
to the presence of a Dzyalozinskii-Moriya chiral term in 
the free energy which impairs the homogeneity of the 
magnetic ordering and gives rise to a continuous helical 
twisting of the ferrimagnetic order parameter with very 
long wavelength. The vanishing of the ME response at 
the low-field phase of Cu20Se03 suggests a multi-domain 
structure of flat helices, where the propagation vectors of 
the different helices are pinned along the preferred axes 
of the system. In the intermediate-field phase the propa- 
gation vectors of the helices are aligned along the applied 
magnetic field H a . For T> 30 K our magnetization mea- 
surements suggest a more complicated state and a possi- 
ble thermal reorientation of the anisotropy. Furthermore 
we also find a narrow range (2K) close to the transition 
temperature T c with a double-peak anomaly in the MEjj 
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APPENDIX 

A. Lifshitz invariants for the MEg effect 

Here we analyze the Lifshitz-type invariants be- 
tween magnetic order parameters and dielectric 
polarization. 4,29 ' 30 These (free) energy terms are linear 
in one spatial derivative of the magnetic order param- 
eter and are allowed in the present compound. As 
we explained above and show also below, these ME 
terms are most likely not related to the observed MEjj 
effect, despite the fact that they give the same angular 
dependence with the term of Eq. (7). However, in a 
systematic expansion of the ME free energy they should 
not be omitted, as they yield terms linear in P and may 
modify certain aspects of the ME behavior, in particular 
in systems with shorter pitch magnetic helices. 
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The Lifshitz invariants are the following seven terms 



B. Magnetocapacitance 



Li = {P,Vi}+ ■ {d,Vj}+, L 2 = {P,»?i}+ • {d,V 3 }- 
L3 = {P,Vi}- ■ {d,Vj}+, ^ = {P,r/J_ • {d,^}. 

L 5 = P x v x d x r 1 x + P y Vid y V V i + P z Vid z Vi 



L 6 = P x rf i d y rf i + P z vtd x tf + P y V fd z ^ 
1.7 = P x riid z Vj + P z ritd y rf + /"'//frV , ( ; 
and their contribution to the free energy is 



8®i 



X n L n 



(32) 



a=l 



where X\ — Xj are the corresponding coupling constants. 
The contribution 5Pl of these terms to the polarization 
can be found again by differentiating with respect to P. 
For example, for E = 0: 



X 



! X if,'' 1 ' '[', + *2T%8»rfi + X^ z d z r, x + X i r 1 z d x r 1 z 
+ X 5V fd x V * + Xt.ifl 0" if + XtrffPrt . 



In order to account properly for the spatial derivatives, 
we may rewrite the position of a given Cu site as = 
r + pi , where p i denotes the relative position of that site 
with respect to the center of the corresponding unit cell 
r. Now, from each term rffWr)^ there will be two types 
of contributions, one proportional to A; 7 and another pro- 
portional to fc 7 (k • Spij). The latter will be disregarded 
in the following since it is quadratic in k. With this in 
mind, and keeping only the uniform component 



one can easily show that only Li and L3 contribute to 
the uniform portion of SP£, and the same happens for 
8P? and 8P£. The final expression for <5Pl is 



8Pl = Xe (A3 - A 2 ) k -ql-qf {e h , e h } 



(33) 



which has exactly the same angular dependence as in 
Eq. (10) above. The difference is that here the pre- 
factor scales with (t]- 1 ) 2 which cannot account for the ob- 
served field dependence, and in particular the sign change 
of P between the "point A" and the 1/2-plateau (see 
Sec. VIII). 



Here we provide a list of invariant terms that are 
quadratic in P and thus may contribute to the MC, in ad- 
dition to the term Uq that we discussed in Sec. IX. These 
invariants are quadratic (in contrast to Uq which is quar- 
tic) in the magnetic order parameters and as such, they 
provide the ~ sin 2 9 angular dependence that is needed 
for the full interpretation of the data shown in Fig. 15(b). 

We restrict ourselves to spatially uniform invariants. 
These are the following five 

U 1 = {r ii ,r )j } + -{P,P} + 

w 2 = K-,^}--{p,p}+ 

% = /';'/;//; + PhH + p 2 v vW 3 

Ui = P 2 x vh V 3 + Ptrgrfi + p^v* 
Us = PlnU) + Phh v 3 + tftfrff 

Including these invariants in Eq. (21) gives 

5 

5$ = $'-$ = -Co Uo - u * ( 34 ) 



where £i — ^6 are the new coupling constants. Repeating 
the steps described in Sec. IX leads to 



d 2 S<I> 
~dPj 



-2Co (vWif ~ 2 (Wry? + „'/;'</f + C 5 ^f) , 



where the second term gives the contributions from the 
invariants U1-U5. Using 



we find that the contribution to the MC from U1-U5 is 



MC' 



[C3(eD 2 + C4(eD 2 + C5(ef i ) 



+ -(C 3 + C4 + C 5 )^r) 2 



±\2 



(35) 



where cry = ±1 is again the sign of rjirjj. 
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